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Abstract 

We consider asymptotically stable scalar-tensor dark energy (DE) models 
for which the equation of state parameter wde tends to zero in the past. The 
viable models are of the phantom type today, however this phantomness is 
milder than in General Relativity if we take into account the varying gravita- 
tional constant when dealing with the SNIa data. We study further the growth 
of matter perturbations and we find a scaling behaviour on large redshifts which 
could provide an important constraint. In particular the growth of matter per- 
turbations on large redshifts in our scalar-tensor models is close to the standard 
behaviour 5 m oc a, while it is substantially different for the best-fit model in 
General Relativity for the same parametrization of the background expansion. 
As for the growth of matter perturbations on small redshifts, we show that 
in these models the parameter 7q = 7' (2 = 0) can take absolute values much 
larger than in models inside General Relativity. Assuming a constant 7 when 
7q is large would lead to a poor fit of the growth function /. This provides 
another characteristic discriminative signature for these models. 
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1 Introduction 



The aim of Dark Energy models is to explain the late-time accelerated expansion of 
the universe [I]. Like for inflationary models producing an early stage of accelerated 
expansion, we have now a wide variety of DE models that can account for the late-time 
background evolution [3J. In the same way that inflationary models are constrained 
by the cosmological perturbations they produce, DE models can be constrained by 
the background evolution and their effect on the growth of perturbations. In principle 
what is basically needed is a smooth component with a sufficiently negative pressure. 
Among DE models, ACDM, although it contains a cosmological constant which can 
be seen as "unnaturally" small, is the simplest model and this model is presently in 
good agreement with observations on large scales (see however e.g. [3]). Another class 
of appealing models are quintessence models containing a minimally coupled scalar 
field [5j. A possible drawback of these models is their inability to violate the weak 
energy condition and to account for a phantom phase. Many more models of ever 
increasing sophistication have been proposed [5]. To make progress, it is important 
to find ways to select those classes of DE models that remain observationnally viable 
and to find characteristic signatures that will enable us to constrain them, or even 
rule them out, with more accurate data at our disposal in the future. 

In particular, an interesting family of DE models are those where gravity is no 
longer described by General Relativity (GR). Indeed there has been considerable inter- 
est recently in DE models with gravitation modified with respect to General Relativ- 
ity, a feature that is quite generic in higher dimensional theories or also in the low en- 
ergy effective action of more fundamental four-dimensional theories. Well-motivated 
models belonging to this class which can be explored thoroughly are scalar-tensor 
dark energy models [TJ El El EDI EH]- Like the usual quintessence models containing a 
minimally coupled scalar field, scalar-tensor models have an additional physical de- 
gree of freedom, namely the scalar partner of the graviton. However, these models are 
more complicated as they have two free fundamental functions in their lagrangian, 
one more function in addition to the scalar field potential. This additional function 
reflects the modification of gravity encoded in the theory. A generic feature of these 
models, like for essentially all alternatives to the cosmological constant model, is that 
DE has a time- varying equation of state. The fact that ACDM fits well a large amount 
of observations should be an incentive to look for models with time-varying equation 
of state still able to compete with ACDM. 

As for all DE models, scalar-tensor DE models are characterized by the accelerated 
expansion they produce at low redshifts but this background effect is common to 
all DE models. The modification of gravity with respect to General Relativity is 
more specific and it expresses itself in particular in the modified growth of linear 
cosmological matter perturbations. In this way, under quite general asumptions, it 
could be possible to determine from the growth of matter perturbations combined 
with the background expansion whether a DE model lies inside General Relativity 
or not [121 [32] (see also [13]). Hence it is very important to investigate how this can 
be extracted from the observations [T5l [16] . Scalar-tensor DE models are certainly 
good examples to investigate this issue. Aside from deeper theoretical motivations, an 
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additional incentive to consider scalar-tensor models could come from the observations 
if these support DE models which have a phantom phase (see also [17] for other models 
that can produce a phantom phase). 

We will consider asymptotically stable, internally consistent, scalar-tensor DE 
models |10| . For such models F — *> F OQ = constant and asymptotic stability is possible 
for wde and fl m — > tt m ,oo- We can consider varying equations of state w EE (z) 
provided wde — ► 0. As we will show these models exhibit a characteristic signature 
in the growth of matter perturbations on large and on small redshifts. This could 
potentially allow us to constrain them or even rule them out. In connection with 
the growth of matter perturbations on small redshifts, these models will illustrate 
results derived in an earlier work [TS], namely the possibility to have a parameter 
7q = -3-(z = 0) whose absolute value is much larger than in General Relativity when 
we write the growth factor as / ~ Q^. Indeed, it was shown in [10J that |7q| < 0.02 in 
models with a constant or a (smoothly) varying equation of state inside GR and hence 
also for ACDM. On the other hand, on large redshifts there can also be an important 
effect on the growth of perturbations. Interestingly, as we will see, we can have models 
for which the growth of matter perturbations on large scales is close to that in ACDM 
and for which the growth of matter perturbations deviates most from that in ACDM 
on small redshifts. Another interesting aspect of a varying gravitational constant is 
its effect on the interpretation of SNIa data and we will show that less phantomness 
is required today by the observations for our models compared to General Relativity. 
Indeed as we consider DE models where wde — >■ in the past, a large amount of 
phantomness is required by the observations if we are inside GR while in our models 
this is much less the case. 



2 An asymptotically stable scalar-tensor model 

It is important to review here the basic aspects and definitions of our scalar-tensor 
models. We consider DE models where gravity is described by the Lagrangian density 
in the Jordan (physical) frame 

L = i (F($) R - Z($) g^d^d^) - U{<$>) + L m (g^) . (1) 

The Brans-Dicke parameter is given by lubd = where a prime denotes a derivative 
with respect to redshift z, and 

G N = (SvrF)- 1 . (2) 

Specializing to a spatially flat universe, the DE energy density and pressure are defined 
as follows 

3F H 2 = p m + p DE (3) 
-2F H = p m + p DE + p DE . (4) 

With these definitions, the usual conservation equation applies: 

Pde = -3H(p DE + p DE ) . (5) 
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With the equation of state parameter wde defined through 

Pde 

w D e = , 

Pde 

the time evolution DE obeys the usual rule 

Pde{z 



= e{z) = exp 



dz 



, 1 + w(z') 
l + z' 



(6) 



(7) 



PDE,0 

Equations ([3]) can be rewritten as 

h 2 (z) = n mfi (l + z) 3 + n DEfi e(z) , (8) 

where Odb,o = 1 — ^m,o by definition as we assume a spatially flat universe. For these 
models 

Pde+Pde = ^ 2 + F-HF + 2(F-F ) H , (9) 

hence the weak energy condition for DE can be violated ( [TDJ , see also [19]). 
We consider viable models satisfying the following requirements: 

1) F — > = constant < F for z — > oo. 

2) The DE equation of state evolves according to wde — * for z — > oo. 

3) Consistency requires /2 = |(^) +|^>0. 

4) We impose ubd,o > 4 x 10 4 . 
Condition 1) is reminiscent of the models considered in [20]. For these models we 



have that 



constant from some redshift on, and different from one 



Gjv,o Gjv,o 

by a few percents only. In order to have accelerated expansion at the present time we 
need some dynamical DE equation of state. We consider a parametrization [21] . [22] 
with a smoothly varying equation of state where DE tends to a scaling behaviour in 
the past 



wde{z) = (-1 + a) + f3 (1 - x) = w + wi 



l + z 



(10) 



where x = — . If we impose wde — > for z — > oo, this prameterization reduces to the 



simple form 



Wde(z) 



l + z 



The corresponding DE evolution reads [21] 

(1 + z) 3 e~ 3 ^ 



elz 



(12) 



For given cosmological parameters fi^^ and fi mj0 the background evolution is com- 
pletely fixed by the parameter (3. We have in particular 



^m. \ Z) 



e 3/3^ ^DE,0 + x 



Clearly, f2 m > ^m,oo (/?) 



^m,oo (/3) 



o 

J 'm,0 



(13) 



(14) 
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The quantity VL m (z) is fixed by fl m ,o an d P- As shown in [TO], the requirement 1) 
implies the inequality (independent of the specific form of F(z)) 

—— > fi m)0 o • (15) 
-TO 

In particular a large amount of phantomness today implies f2 m)00 close to one and 
hence F^ close to Fq. We use further the following ansatz 

This ansatz satisfies exactly F% = |r(z = 0) = while F 2 = ^r(z = 0) = —10(1 — 
^) < 0. For the ansatz ffT6]) . we must have in addition < l+^f2 OE)0 (l— /9) which 
comes from the requirement 0' 2 > at z = 0. So we impose the functions h(z) and 
-j^(z) from which all other quantities can be reconstructed and we check the physical 
consistency for each reconstructed scalar-tensor model. As h(z) is imposed, so is the 
background dynamics. In this way we can compare different DE models inside and 
outside General relativity possessing the same background evolution. Clearly, as the 
DE equation of state parameter w^e tends to zero as z increases, it must start being 
phantom today if it is to pass the observational constraints. This can be seen more 
quantitatively using the constraint on the shift parameter. We see from Figured] that 
the viable models are of the phantom type today (/3 > 1). The possibility to have a 
phantom DE sector today is actually an attractive feature of scalar-tensor DE models 
and is not excluded by the observations. 



3 Some observational constraints 

To get some insight into the parameter window for viable models, we constrain them 
using Supernovae data, BAO (Baryonic acoustic oscillations) data and CMB data. 
We have to maximize the probability function 

P(fi m , ,/3)oce-^ 2 . (17) 

where \ 2 — Xsn + Xa + Xr f° r a background evolving according to eqs.([51). f[T2l . We 
use a sample consisting of 192 Supernovae [6l [2] for which 



192 



with 



i2 



, ( , \ f Zi dz\ 15 , G cS (zi) , . 

m h>i = 5 log ^(1 + Zi ) ^ — J + no + — log , (19) 

where /i = 25 + 5 log ^ c ^° e j , the distance modulus \i is the difference between the 

apparent magnitude m and the absolute magnitude M. The important quantity G c q 
is defined as 

F + 2(dF/d^_ 1 + 2u- B l D 

GeS ~ f + \{dF/d*)* ~ Gn TTf^ ' { ) 
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We have for ui BD ^> 1 



G 



cff 



G 



N 



[1 + -U 



1 

— L 
2 



bd) 



(21) 



In particular G c q$ ~ Gn,q due to the well-known strong solar system gravitational 
constraint ubd,o > 4 x 10 4 . We get rid of the nuisance parameter H using the simple 
way suggested by [23], integrating over H gives essentially the same result. Note the 
addition of the last term in eq. (1191) which takes into account a varying gravitational 
constant [23]. This term allows to discriminate different scalar-tensor models using 
SNIa data even for similar background expansion. We will come back to this point 
below. 

The BAO constraints can be expressed as a constraint on the quantity A 



A{z) 



h(z) 



with [25] 
and 



A = 0.469 ±0.017 



2 _ (A{z = 0.35, fi m , , /?) -0.469) : 
Xa — 



:22) 



(23) 



(0.017) 2 

We have finally a constraint on the shift parameter extracted from the CMB data 

,1089 dz 



(24) 



/•i 
Jo 



with [26] 
and 



h{z) 

R= 1.70 ±0.03 , 



2 (R(n mfi ,f3)-1.7f 
X R = ~ 



(0.03)' 



(25) 

(26) 
(27) 



As we can see from Figure 1, the shift parameter constrains our model to be of the 
phantom type today, (3 > 1. This is expected because in our model the equation of 
state of DE tends to that of dust in the past. If we remember that a cosmological 
constant agrees fairly with the data, our model must compensate by being phantom 
on small redshifts. 

We would like now to look more closely at the effect of a varying (effective) grav- 
itational constant on the measurement of luminosity distances. Let us write eq. ([T9l 
in the following way 



Hth,i = [1 + Q(Zi)\ 5 lo S ( (1 ± Zi 
where we have introduced the quantity 



dz 
~h 



(2f 



i G c (r(z) 

3 log J tt v > 

4 log [(1 + tt ¥1 ' 



(29) 
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Foo/ F 




P 






c 


Pi 




7o 


7o 


0.93 


n qi+u.u4 

u - oi -0.04 


1 09 +u ' iy 

i - UJ -0.13 




Q2 +u - 04 

u - JZ -0.04 


n qq+0.04 

u - yj -0.04 


u - yy -0.03 




u - d ^-o.oi 


-0.07™ 


0.94 


r, on+0.04 


i -i O+0.12 
i - iz -0.13 




0.93I™ 4 


0.98^ 


u - yy -0.02 




^4+ - 01 
u - d ^-o.oi 


-0 06+ - 03 

u - uu -0.02 


0.95 


-?o+ 04 


i 15+0-13 




0.93t^ 4 


o.98ig : gS 


u - yy -0.02 




u - oq: -0.01 


-n n4 +0 - 02 

u - u ^-0.01 


0.96 


29+ - 04 

u - zy -0.03 


1.181°;- 




U - JO -0.04 


q7+ 04 

u - y ' -0.04 


q8+°- 02 




o ^ +0 - 02 

u - dd -0.01 


-0 m+ - 05 


GR 


n 97+0.04 
u - z ' -0.03 


1.281°;^ 




q^+ 04 

u - Jd -0.03 


u - yj -0.03 


n q7+ - 02 

U - J ' -0.02 




o % +0 - 01 


01+ 002 

u - u± -0.002 



Table 1: We summarize in this table the best-fit models for given parameter when 
all data are taken into account, with 2a errors. The last line corresponds to General 
Relativity (a constant gravitational constant). It is seen that a varying gravitational 
constant G c s can have a nonnegligible effect. In particular, though the quantity 
f2 m>00 is higher in GR than in the scalar-tensor models, due to the value of F^, C 
is closer to 1 in these models. Hence the growth of matter perturbations on large 
redshifts for these models, 5 m oc a Pl , is closer to the standard one as in ACDM, see 
eqs.([37D,(l38D,(j39]). 



We can write in full generality 

G eS (z) 



G 



cS,0 



A(z) . (30) 



In the models studied here, A(z) is positive definite and at most of the order of 
a few percents, A < 0.07, while A(z = 0) vanishes by definition. For any viable 
scalar-tensor model we have on very small redshifts 



1 . , , ,-i v , fFo 



A(z « 0) ~ - {ugh-^BDja) + [-J ~ lj « • (31) 

We can give a more formal expression for A(z) but for our purposes it is not needed 
here. Whenever the quantity A is small, we obtain 

^ log ~ 1.63 A A < 1 . (32) 

and hence also 

|£(2«0)| -0.33 - r «l. (33) 

I logz| 

Hence it is seen from ( 1281 . (j3Tj) . (1331) that the effect of a varying G efT for scalar-tensor 
dark energy models is negligible on very small redshifts. We have indeed checked it 
with our models using SNIa data on redshifts z < 0.05. Furthermore, it is also clear 
from (13"2"1) that this effect cannot be very large whenever A <C 1. 

Nevertheless, as we can see from Table 1, a varying gravitational constant G e g, 
characterized in our models by the parameter F^, can have a nonnegligible effect. In 
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Figure 1: a) On the left, the shift parameter R (see eq. (125jl ) is shown in function of 
(3 for several values of Q m ,o- The interval of viable models all correspond to phan- 
tom DE, (3 > 1 which was expected because wde —> in the past, b) On the 
right, the quantity — (with arbitrary normalization, i.e. — | z =ioo — 1) is shown 
for various cases: General Relativity (GR) (blue), ACDM (red) and scalar-tensor 
models (green) characterized by the value of F^. We have from top to bottom 
= 0.9301, 0.94, 0.95, 0.96, 0.97, 0.98. The first value corresponds toCwl and the 
linear growth on large redshifts is essentially similar to ACDM. Except for ACDM, 
all the models displayed here have the same background behaviour with Q m ,o — 0.3, 
and Q m ,oo = 0.93. 



particular it is seen that our models with < F can fit the data with dark energy 
which is less of the phantom type today than it would have to be in the corresponding 
DE model inside GR with a background expansion of the type ([HD, (fT2"j) . 

4 Linear growth of perturbations 

Let us turn now to the dynamics of the linear matter perturbations. As shown in [8], 
these perturbations satisfy a modified equation of the type 

5 + 2H5 - AttG cS p m S = , (34) 

with G e s given by eq. (1201) . Equation ( 1341) can be seen as a minimal modification to 
the growth of linear perturbations which comes from the modification of Poisson's 
equation 

^ <j> = -4tt G p 5 -> \ <j) = -4tt G cS p 5 . (35) 

This modification reflects the fact that the effective coupling constant describing the 
gravitational interaction of two close test masses is given by G c g. This is so on all 
cosmic scales of interest where the dilaton field is essentially massless. It should be 
stressed that the modification in ( 1341) is scale-independent and can appear in many 
modified gravity models (see e.g. [27]), it can even appear in DE models inside 
General Relativity if one is willing to consider DE with unusual properties [28J. 
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It is convenient to introduce the quantity / = 4jj^, the growth factor of the 
perturbations. In function of /, the linear perturbations obey the equation 

df 2 1 f d\nQ m \ 3 G eS 



^♦H 1 -^'-^*"- (36) 

with x = In a. The quantity 5 is recovered from /, 5(a) = 5i exp f(x')dx' . We 
see that / = p when 5 oc a p , in particular / — > 1 in ACDM for large 2 while / = 1 
in an Einstein-de Sitter universe. In our model, G e ff — ► Ges t00 , — ► ^m,oo, these 
quantities tend rather quickly to their asymptotic value for z ^> 1. 
Introducing the quantity C with 

< C = — &<m,oo = ~^~^m,oo < 1 > (37) 

we see that in the asymptotic regime G e fr — > G e fr,oo = -f^-, the perturbations obey a 
scaling behaviour 

5 = D^ 1 + £> 2 a P2 , (38) 



with 



Pi = ^(-1 + V1 + 24C) (39) 



p 2 = -(-1 - Vl + 24C) . (40) 

As pi < 1, we see that the growing mode of the perburbations allways grows slowlier 
than in a ACDM universe (or in an Einstein-de Sitter universe) for z ^> 1. Therefore 
the amplitude of the linear matter perturbations on small redshifts before formation of 
structure starts, compared to the amplitude of perturbations derived from the CMB 
data, can be significantly different, and smaller, from that in ACDM. In particular, 
the perturbations will grow nonlinear on lower redshifts, structure formation starts 
later. Further, the bias h derived from o"§ should be larger in these models than it 
is in ACDM. On the other hand, for a model for which C is very close to 1 (but 
still satisfying C < 1), both linear perturbations modes will evolve essentially like 
in a ACDM universe untill low redshifts where a significant departure can appear. 
Another important issue is to characterize this departure on small redshifts. 

It is well known that for a ACDM universe it is possible to write / ~ where 
7 is assumed to be constant, an approach pioneered in the literature some time ago 
[291 EH]. There has been renewed interest lately in this approach as the growth of 
matter perturbations could be a decisive way to discriminate between models that 
are either inside or outside General Relativity (GR). Clearly it is possible to write 
allways 

/ = ■ (41) 

For ACDM we have 7 = j(z = 0) w 0.55. As was shown in [18], for ACDM we 
have 7g = ^| « -0.015. For £l mfi = 0.3, 7 = 0.555, slightly higher than the 
constant ^ = 0.5454 derived in [31] for a slowly varying DE equation of state and 
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Q m ~ 1. There is a very slight difference on small resdshifts z < 0.5 between the 

_6_ 

true function /a(<z) and Qm , one could as well use 7 = 0.56 and the agreement would 
be even better. As Q m — > 1 the differences are important only on small redshifts. 
Note also that we find a slightly negative slope 7 so that 7 comes closer to as 
z increases. A definite departure from these values could signal a departure from a 
ACDM universe. More importantly as we will see later a large value for 7 could be 
a hint for a DE model outside GR. This aspect was already emphasized in [18j. We 
illustrate in Figure 2 the behaviour of 70 and 7 in function of ft and Q m ,o inside GR. 
As emphasized in [18], equation (1361) yields the following identity 



7q 



-n£o - 3(7o - \) w cS , + ~ niy - i 



(42) 



whenever g e ^'° = 1 to very high accuracy, which is certainly the case in scalar-tensor 
models as we have seen in Section 2. In other words we have a constraint 

/(7o, 7o> fi m,o, w de ,q) = . (43) 

As was shown in [18], this constraint takes the following form 

7 ~ -0.19 + d (70-O.5) d^3. (44) 

The coefficient d depends on the background parameters d = d(wDE,0i^m,a)- F° r 
given background parameters fl m ,o and wde,o, 7 will take the corresponding value 
7 (7o). The value of 70 realized will depend on the particular model under consider- 
ation and can be obtained numerically. Typically we will have 7 7^ 0. For models 
inside General Relativity |7 | < 0.02 was obtained. For example for constant wde, 
7 is allmost independent of wde = wde.Oi with 7 rs —0.02 for O m>0 = 0.3, while 
at the same time 70 can have a nonnegligible variation, we have 0.545 < 70 < 0.565 
for —1.4 < WoEfl < —0.8 and f2 m0 = 0.3. We can compare scalar-tensor models 
with different values of but identical background evolution parametrized using 
the parameters ft and fl m ,o, according to eqs.([8]), ffTTl) . f|T2l . These models can be 
distinguished in all observations affected by a varying gravitational constant but they 
will be undistinguishable with respect to purely background constraints. In partic- 
ular, they yield different perturbations growth factor /, or equivalently different 7. 



All models shown in Figure [3] have exactly the same background expansion so that 
the difference in the growth of matter perturbations is solely due to the modification 
of gravity. As we can see from Figure [31 some models can be easily distinguished 
from each other using the growth of matter perturbations through the difference in 
the behaviour of 7(2) on small redshifts. We note also for the models displayed 
in Figure [3] that the model inside GR is easily distinguished from the scalar-tensor 
models having a large slope 7 . For these ST models we have that 7 is negative and 
it can be large, while it is (slightly) positive for the same background evolution inside 
GR. Actually as we can see from Figure 2 7 > is true inside GR for the interesting 
range of cosmological parameters ft and f2 mi o, while we have generically 7 < for 
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Figure 2: a) On the left, 70 is shown in function of the parameters £l m ,o and /3 for 
the ansatz ffTTj) assuming General Relativity, b) On the right, the corresponding 7q 
is displayed and it is seen that \j' \ < 0.015. 



1.1 




0.5 1 1.5 2 0.5 1 1.5 2 



z z 

Figure 3: a) On the left, the function f(z) is shown corresponding to the models of 
the right panel of figure [lj The growth factor is constrained by observations of the 
galaxy bias factor b and the redshift distortion parameter j3 [32] via / = b (3. The 
data are given here for reference. We see that one needs to go at redshifts z > 0.5 in 
order to be able to distinguish between the different models. It is also clear that the 
present observations are not discrimative. Though not obvious from the behaviour of 
f(z), the behaviour on small redshifts is very different as seen on the right panel, b) 
On the right, the function 7(2) is displayed corresponding to the same models as on 
the left panel. But the order of the scalar-tensor models is reversed here, we have now 
from bottom to top = 0.9301, 0.94, 0.95, 0.96, 0.97, 0.98. It is seen that 7O) has 
a quasi-linear behaviour on small redshifts and that large slopes can be obtained for 
scalar-tensor models when C is close to 1. Interestingly, this is also when the linear 
growth of perturbations is essentially the same as in ACDM for z 3> 1 in the matter 
dominated stage. 
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Figure 4: In green the growth function f(z) is shown corresponding to the scalar- 
tensor model with ^ = 0.9301 while the background is parametrised by fi m;00 = 
0.93 and f2 mj o = 0.3. This function is compared to three fits with a constant 
7 : 0.4, 0.48, 7 = 0.53, where 7 is the value at z = of the true j(z). For 
7o = 0.53, the fit is very good (by definition) on very small redshifts but bad on 
larger redshifts. Choosing 7 = 0.4 gives a good fit on large redshifts but is bad on 
small redshifts. Taking 7 = 0.48 yields a fit that is clearly different from /. 



the interesting scalar-tensor models. Also, as already noted in [J8], |7g| < 0.02 inside 
GR while we see that it can be substantially larger outside GR. Obviously, when |7q| 
is large, assuming 7 = constant yields a poor approximation. Another interesting 
feature is the quasi-linear behaviour of 7(2) on small redshifts < z < 0.5. Such a 
behaviour could be probed observationnally and could allow to discriminate models 
whose parameter 70 are close to each other. The lower the value of 70, the better this 
potential resolution which improves as well when wde,o decreases. 

This picture remains essentially the same for the best-fit models of Table 1. It 
is seen that the growth of matter perturbations on large redshifts is closer to the 
standard behaviour 5 m oc a (C = 1) in the best-fit scalar-tensor models of Table 
1 than it is for the model inside GR with same parametrization of the background 
expansion. For the latter model, though a very pronounced phantom behaviour is 
needed today (see Table 1), it is not sufficient to make fi mj0 o higher than 0.95 which 
has a crucial effect on the growth of matter perturbations on large redshifts. 

To summarize, we have shown that our models have a characteristic signature 
in the growth of linear matter perturbations. On large redshifts inside the matter- 
dominated stage, we find a scaling behaviour for the matter perturbations which can 
substantially differ from ACDM and also from General Relativity (GR) with identical 
background evolution characterized, besides Q m ,o, by fi mi00 or (3. On small redshifts 
we find again a possible significant departure from ACDM and models in General 
Relativity (GR) with same parametrization of the background expansion. Even for 
those models in which the growth of matter perturbations on large redshifts is close 
to that in ACDM, we find a large (negative) slope 7^ with |7q| much larger than 
in GR, whether ACDM or GR with an identical background evolution. For these 
models assuming a constant 7 would necessarily lead to a poor fit of the growth 
function / (see Figure 4). Such a behaviour on small redshifts would constitute a 
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characteristic signature of our DE model being outside GR. Interestingly, those models 
that mimic ACDM on large redshifts are most easily distinguished from ACDM on 
small redshifts through their slope 7g. Though the results derived here are to some 
extent model dependent, it is clear that the growth of matter perturbations, especially 
when combined on small and large redshifts, can efficiently probe the nature of Dark 
Energy and in particular help in assessing whether we are dealing with a modified 
gravity DE model or not. 
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